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Explicit Frames for Deterministic Phase Retrieval via PhaseLift 

Michael Kech^’* 

^Department of Mathematics, Technische Universitat Munchen, 8574-8 Garching, Germany 

(Dated: January 24, 2017) 

We explicitly give a frame of cardinality 5n — 6 such that every signal in C” can be recov¬ 
ered up to a phase from its associated intensity measurements via the PhaseLift approach. 
Furthermore, we give explicit linear measurements with 4r(n — r) -|- n — 2r outcomes that 
enable the recovery of every positive semidefinite n x n matrix of rank at most r. 
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I. INTRODUCTION AND MAIN RESULT 

Phase Retrieval is the task of reconstructing a signal x G C"" up to a phase from intensity 
measurements. 

In [1] it was shown that m > 4n — 2 generic intensity measurements suffice to discrimi¬ 
nate any two signals in C”" up to a phase. With a similar approach this result was slightly 
improved to m > 4n — 4 in [2] ^. The bound m > 4n — 4 is known to be close to optimal. 
More precisely, by relating phase retrieval to the problem of embedding complex projec¬ 
tive space in Euclidean space, it was shown in [6] that, up to terms at most logarithmic 
in n, m > 4n — 4 intensity measurements are necessary to discriminate any two signals in 


‘Electronic address: kech@ma.tum.de 

^ In the context of pure state tomography, [3-5] show that the 4n — 4 bound also holds for von Neu¬ 
mann measurements. In addition similar bounds for the recovery of low-rank matrices with constrained 
measurements are provided in [3]. 
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C"" up to a phase. However, [1, 2] do not provide a tractable recovery scheme. A result 
indicating that some redundancy is needed in order to allow for computationally efficient 
phase retrieval is given in [7]. 

There have been several approaches that do provide recovery schemes [8-10], in the 
present paper however we focus on the approach of [11] known as PhaseLift. Their approach 
consists of two steps: First, phase retrieval is lifted to the problem of recovering rank one 
Hermitian matrices from linear measurements. Secondly, by means of a convex relaxation, 
the recovery problem is formulated as a trace norm minimization over a spectrahedron. 
The authors of [11] then prove that 0(n) intensity measurements suffice to recover a 
signal modulo phase with high probability by solving the relaxed optimization problem. 
Furthermore, stability guarantees for the recovery were established in [12, 13]. While 
these convex relaxations are in principal tractable, solving them becomes computationally 
expensive with increasing signal dimension [14]. 

However, [11-13] still leave room for improvement. For example, by working with 
Gaussian random vectors additional structure that might facilitate the use of PhaseLift is 
not incorporated and also from a practical point of view Gaussian random vectors might 
not be desirable. Recently, it was shown that a partial derandomization of PhaseLift 
can be achieved by using spherical designs [15, 16]. The purpose of the present paper is 
similar. However, rather than drawing the measurements from a smaller, possibly better 
structured set, we aim for finding explicit measurements that allow for phase retrieval via 
PhaseLift. Another deterministic approach to the phase retrieval problem was introduced 
in [17]. They improved their results in [18], providing recovery algorithms together with 
explicit error bounds for phase retrieval with 6n — 3 frame vectors. 

Our contribution is the following: We explicitly give 5n — 6 intensity measurements 
from which every signal in can be reconstructed up to a phase using PhaseLift. More 
precisely, for A: G {1,..., 2n — 3} let 

Vk := (^1, Xk e^, xl ..., , Xfc G M \ {0}. (1) 

Furthermore denote by {ej}j£{o,...,n-i} the standard orthonormal basis of C”. 

Theorem I.l. // xi < X 2 < • • • < X 2 n- 3 , then every signal x G C” ean he reeonstrueted 
up to a phase from the 5n — 6 intensities 

{|(eo,x)|^,...,|(en_i,x)p,|(ni,x)p,|(l;i,x)p,...,|(n2n-3,a:)P,|(F2n-3,a:^)|^} 


via PhaseLift. 

This result is stated more carefully in Section HI as Gorollary HI.2. Its proof relies on 
the results of [19]. 

Let us highlight three features of this result: 

1. To our knowledge the 5n — 6 is the smallest number of intensity measurements that 
allow for a uniform and computationally tractable recovery. 

2. Results based on random intensity measurements typically guarantee that the recov¬ 
ery succeeds with high probability if the number of measurements exceeds a given 
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threshold which is usually determined up to a multiplicative constant. As opposed 
to this, Theorem I.l comes with two advantages that might be desirable from a prac¬ 
tical point of view: First, the recovery is not just guaranteed to succeed with high 
probability but indeed works deterministically. Secondly, since the measurements 
are given explicitly there is no need for finding a suitable value for the threshold. 

3. Theorem I.l merely requires 5n — 6 intensity measurements. This illustrates that n 
additional measurements as compared to the nearly optimal bound of [1] suffice to 
render PhaseLift feasible. 

The approach we take originates from low-rank matrix recovery [20-24] and indeed 
the previous results can be generalised to this setting: In Section III, we give an explicit 
family of linear measurements with 4r(n — r) -|- n — 2r outcomes from which every positive 
n X n matrix of rank at most r can be recovered by means of a semidefinite program. 
This strongly relies on the construction of the null spaces of such measurements given in 
[19]. Our contribution is to explicitly characterize the orthogonal complements of these 
null spaces leading to the proofs of our main results. 

Finally we also prove a weak stability result in Section IV, showing that the reconstruc¬ 
tion error is linear in the error scale. As we do not know how to estimate the constant of 
proportionality appearing in the stability bound, this result is not of practical relevance, 
but might give a roadmap for proving stability in the future. However, we provide some 
numerical results that might indicate the constant’s qualitative behaviour. 


II. PRELIMINARIES 

Let us first fix some notation. By M{n,q) {M{n,q,M.)) we denote the set of complex 
(real) nx q matrices. The transpose (conjugate transpose) of a matrix A G M{n,q) is 
denoted by (A*). For i G {0,... , n — 1}, j G {0,..., (7 — 1}, we denote the entry in 
the i-th row and j-th column of a matrix A G M{n,q) by Aij By H{n) we denote the 
real vector space of Hermitian n x n matrices. We equip H[n) with the Hilbert-Schmidt 
inner product and || • ||2 denotes the Frobenius norm. By 5" we denote the set of positive 
semidefinite n x n matrices and by C we denote the subset of positive semidefinite 
matrices of rank at most r. In the following we assume that rG{l,...,[n/2]—1}^. The 
set of linear maps M : H (n) —)■ M™' is denoted by A4 (m). 

Definition II. 1. (m-measurement.) An m-measurement is an element of 

In the following we denote an m-measurement simply by measurement if we do not 
want to specify m. 

Remark For each m-measurement M there exists a unique G := {Gi,... ,Gm) G H{n)^ 
such that 

M{X) = (tr(GiA),...,tr(G^A)) 


Note that the indices we use to label matrices begin with 0, not with 1. 
^ [fc] := the smallest integer i such that i > k. 
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for all X G H{n). By Mq we denote the m-measurement associated in this way to an 
G G H{n)'^. In the following we sometimes use this identification to speak of elements 
G G H{n)^ as m-measurements. 

Definition II.2. (r-complete.) A measurement M is ealled r-eomplete iff M{X) M{X') 
for all X G Sf and X' G with X 7 ^ X'. A tuple G G H{n)"^ is ealled r-eomplete iff Mq 
is r-eomplete. 

Given a measurement M and a measurement outcome b = M{X), X G Sf, consider the 
following well-known semi-definite program [20, 22, 23] ^ 

minimize tr(y) 

(2) 

subject to y > 0, M(y) = b. 

The significance of the r-complete property is due to the following observation: 

Proposition II. 1. Let M be an r-eomplete measurement and let X G Sf. If b = M{X), 
then X is the unique minimizer of the semidefinite program (2). 

Proof. Let X G Sr(C^) be a Hermitian matrix of rank at most r and let M be an r- 
complete measurement. Then, X is the unique feasible point of the spectrahedron 

{y G II(n) : y > 0, M(V) = M(X)}. (3) 

This follows immediately from {y G iL(n) : y > 0, M{Y) = M{X)} = {Y G 5” : 
M(Y) = M(X)} and the definition of r-complete. □ 

Remark Note that if 1 G Range(M*), the r-complete property also is necessary for a 
deterministic reconstruction via the semidefinite program ( 2 ). 

This shows that for an r-complete measurement the semidefinite program (2) reduces to a 
feasibility problem. 

Finally, let us state the observation of [19, 25] which gives a useful characterization of 
the r-complete property: 

Proposition II.2. A measurement M is r-eomplete if and only if every nonzero X G 
Ker{M) has at least r -|- 1 positive eigenvalues. 

Proof. Consider the set A:={y — y:y G Sf, Z G S^} and note that every X G A has 
at most r positive eigenvalues. Furthermore, note that a measurement M is r-complete if 
and only if A n Ker(M) \ {0} = 0. 

Now assume that every X G Ker(M) \ {0} has at least r -|- 1 positive eigenvalues. 
Since every Y G A has at most r positive eigenvalues we find Y 4 Ker(M) \ {0|, i.e. 
AnKer(M)\{0} = 0. 

Conversely, assume that M is r-complete. A clearly contains all matrices with at most 
r positive eigenvalues and hence Ker(M) \ {0} cannot contain an element with r or less 
positive eigenvalues. 

□ 

This is a convex relaxation of the rank minimization problem. 
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Remark If every nonzero X G Ker(M) has at least r + 1 positive eigenvalues, then every 
nonzero X G Ker(M) also has at least r + 1 negative eigenvalues since X G Ker(M) implies 
-X G Ker(M). 


III. RECONSTRUCTION OF LOW-RANK POSITIVE MATRICES 

Our approach relies on [19] where a method to construct the null spaces of r-complete 
m-measurements for m = ^r{n — r) + n — 2r \s provided. Their construction is based on 
the ideas of [26], details can be found in Appendix A of [19]. 

First, we focus on the phase retrieval problem. 

Theorem III.l. Let 

vivl Viv\ V 2 n-:iV* 2 n-Z V2n-3V2n-3 \ 

\\v2n-3V2n-3\\2' ll^2n-3W2n-3 lb / ’ 

where the Vi are defined in Equation ( 1 ). If xi < X2 <■■■ < X2n-3, then the measurement 
Mg is 1-eomplete. 

The proof of this theorem can be found in Section V. 

Remark From the proof of this result it is easily seen that the kernel of Mq is independent 
of the choice of the Xi. Thus, for the purpose of robustness, the Xi should be chosen such 
that the smallest singular value of Mq is maximized. 

Let us next state Theorem I.l in a more precise way. 

Corollary III. 2. (Phase Retrieval via PhaseLift.) Let M be a measurement given by 
Theorem III.l and let x G C”. If b = M{xx*), then xx* is the unique minimizer of the 
semidefinite program (2). 

By Proposition II. 1, this is an immediate consequence of Theorem III.l. 

Let us next focus on the recovery of low-rank positive matrices. This, however, requires 
some further definitions: First, let 

Cf := {X G H{n) : triXae*) =0, 2r - 1 < i + j < 2(n - r) - 1, i / j, }. (4) 

E.g. C'f C II{n) is the subspace of n x n diagonal matrices and Cj is the subspace of iL(7) 
of the from 


G . I 6q6q, • • • 5 1^72—1’ 


/*****oo\ 

* * * * 0 0 0 

> 1 = * * 0 0 0 * 

* * 0 * 0 * * 

* 0 0 0 * * * 

0 0 0 * * * * 

\0 0 * * * * */ 
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For X G M\{0}, k G {2r— 1,..., 2(n —r) — 1}, define the Hermitian matrices Rk{x), Ik{x) G 
{C^)^ by ^ 

{Rk{x))ji := 5j+i^kx\ j,l £ {0,... ,n- 1}, j > I, 

{Ik{x))ji := i5j+i^kx\ j,l £ {0,... ,n- 1}, j > I, 

where 5ij denotes the Kronecker delta. E.g. for n = 5, r = 2 these are 

/OOOION /OOOiON /0 000 1\ /0 000 i\ 

OOa^OOj 0 0 ix 0 o] (OOOa^Oj 0 0 0 ix 0 \ 

R 3 (x) = 0 2 ! 0 0 0 , 73 ( 2 ;)= 0 -ix 0 0 0 ,H 4 ( 2 !)= 0 0 0 0 0 , 74 ( 2 :)= 0 0 0 0 0 , 

10000 -iOOOO 02:000 0 -ix 0 0 0 

\0 0 0 0 0 / VOODOO/ \1 0 0 0 0 / \-i 0 0 0 0 / 
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Theorem III. 3. Let Go be a basis of Cf and let xi,X 2 , ■■■ ,Xr G M \ {0} with xi < X 2 < 
... < Xr- For k £ {2r — 1,, 2(n — r) — 1} define 

Gk {Ik{xi)j Rk{xi)j • • • ) Lk{Xr), Rk{Xr)) ■ 

and let G := GqU G 2 r-i U • • • U G 2 (n-r)-i Then the measurement Mq is r-eomplete and 
|G| = 4r(n — r) + n — 2r. 

Remark If an m-measurement is injective when restricted to Sf, it was shown in [6, 27] 
that, up to terms at most logarithmic in n, we have m > 4r(n — r). Furthermore, in [3, 6] 
it was shown that there indeed exist injective m-measurements for m = 4r(n —r). Thus, it 
might be worth noting that the measurements given by Theorem III.3 solely require n — 2r 
additional measurement outcomes as compared to the nearly optimal bound 4r(n — r). 

Finally, by Proposition II.1, the measurements given by Theorem III.3 allow for the recov¬ 
ery of low-rank positive matrices. 

Corollary III.4. (Recovery of low-rank positive matrices.) Let M be a measurement given 
by Theorem III.3 and let X £ Sf. If b = M{X), then X is the unique minimizer of the 
semidefinite program (2). 

IV. STABILITY 

In this section we discuss the stability of r-complete measurements. 

Assume there is an error term E £ H(n) that perturbs the matrix Xr £ 5” we intend to 
recover to the matrix X = Xr -£ E. Measuring with an r-complete measurement M yields 

® As Rk{x), Ik{x) G (C")^ both have vanishing diagonal and since they are hermitian, it suffices to define 
all elements above the diagonal. 

® For tuples of Hermitian matrices X ~ {X\, ..., Xi) G Hfii)', Y ~ (Yi,..., Yj) G H{n)^ we define their 
union X U Y to be the tuple X U Y := (Xi ,..., Xi,Yi, ■ ■ ■ Xj) £ H{n)'+L 
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the perturbed outcome b = M{X). Clearly, the matrix Xr cannot always be perfectly 
recovered from the outcome 6, however, if ||M(ill )||2 is small, there is a recovery procedure 
that yields a matrix close to Xr- For that purpose, consider the following well-known 
optimization problem 

minimize tr(y) 

(5 

subject to y > 0, ||M(y) — b \\2 < e 
where e > 0 is a constant representing the error scale. 

Theorem IV.1. (Stable recovery of low-rank positive matrices.) Let M be an r-complete 
measurement and let e > 0. There is a constant Cm > 0 independent of e such that for all 
Xr £ Sf and E £ H{n) with ||M(F1)||2 < e, any minimizer Y of (5) for b = M{Xr + E) 
satisfies 


||y — Xr \\2 < Cm^- 

Remark In the proof of this theorem we show that Cm < + “) where amin is the 

•^min ^ ' 

smallest singular value of M and k := —vaceKz^Ker{M),\\z\\ 2 =i^ri-r{Z) However we do 
not know how to compute k for a given r-complete measurement M and hence we cannot 
make this bound more explicit. 

The proof of this theorem can be found in Section V. In order for this result to be a practical 
stability guarantee, one would have to estimate the constant Cm- At this point we do not 
know how this can be achieved. In order to indicate the magnitude of the constant Cm, 
let us next present some numerical results. For this purpose consider the tuple 

7l(l) -Rl(l) 7?2n-3(l) \ 

’ ||/2n-3(l)||2 ’ ||i?2n-3(l)||2y' 

and note that by Theorem III.3 the associated measurement Mg„ is 1-complete. Figure 
1 presents numerical results that might indicate the scaling of Cmg„ the sequence of 
measurements {MG„)n&n- 

Just like in [12], this recovery scheme can also be used for the phase retrieval problem. 
For a Hermitian matrix A £ H{n), we denote by Eig(yl) G the tuple of eigenvalues of 
A ordered decreasingly together with their multiplicities. Furthermore, we define Aj(H) ; = 
Eig(H)i_i, i G {l,...,n}. 

Proposition IV.2. (Stability for Phase Retrieval.) Let X = xx* + E, where x £ C"' is 
the signal and E £ H{n) is an error term. Let M be a 1-complete measurement and let 
e y ||Af(Fl)|| 2 . Furthermore, let Y be any minimizer of the optimization problem (5) for 
b = M{X) and set x := y^~X(jY)x' where x' £ S'^~^ is an eigenvector ofY with eigenvalue 
Ai(y). Then 


Gn := 


eoeo, 


I ^n—l^n— 


XX* — XX*||2 < 2Cm^, 


Xn-r is defined later this section. 
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FIG. 1: For each n G {3,6,..., 75} we choose uniformly at random a normalized vector x G C" 
and an error term / G with ||/||2 < e := 10“^. Then we run the program (5) with the 

outcome b = Mc„{xx*) + /• The figure shows the maximum value of ||F — xx*\ilt for 2200 
repetitions where Y is the minimizer of (5). 


where Cm is the constant given by Theorem IV. 1. Furthermore, for some ip G [0, 27r) we 
have 


2V2Cm 



This result follows from Theorem IV. 1 by a straightforward computation. The proof is 
given in Section V. 


Remark The proofs of V.5 shows that the above stability results also hold true the fol¬ 
lowing recovery scheme: 

minimize ||M(y) — b \\2 

( 6 ) 

subject to y > 0, 


where M is r-complete and b = M{Xr + E), Xr G 5”. 


V. TECHNICAL APPENDIX 

Let us first introduce some notation we use throughout this section. Let A G M{n,q), 
i G {0,...,n — 1}, j G {0, ... ,q — 1}. By A^i we denote the (n — 1) x q matrix obtained 
from A by deleting the i-th row and by A’^ we denote the n x (g — 1) matrix obtained 
from A by deleting the j-th column. By A{i} we denote the i-th row of A and by A[j] we 
denote the j-th column of A. Furthermore, for k G {0, ... ,n + q — 2}, we denote the k-th. 
anti-diagonal of A by A{k), i.e. A{k) := 


The ordering is such that the matrix element with smaller i comes first. 
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A. Proof of Theorem III.3 

Since Theorem III.l is obtained by manipulating the measurements obtained from The¬ 
orem III.3 we begin by proving the latter. The construction we give in the following yields 
a more general class of r-complete measurements than the ones given by Theorem III.3 
and it strongly relies on the notion of totally non-singular matrices. 

Definition V.l. (Totally non-singular.) A matrix A G M{n,q) is called totally non¬ 
singular if A has no vanishing minor. 

The following lemma is a central ingredient for the construction given in the following. 

Lemma V.l. Let q G {l,...,n — 1} and let A G M{n,q) be totally non-singular. Then, 
there exists a totally non-singular matrix B G M{n,n — q) sueh that A*B = 0. 

Proof. We give a proof by induction in the dimension n for q fixed. 

Base case. Let us begin with the base case n = q -\- 1. Note that for a given A G 
M(q -\-1, q) there always exists a nonzero matrix (actually just a vector) B G M(q -\- 1,1) 
such that A*B = 0, in particular if A is totally non-singular. 

Since B exists, it is enough to prove that if A is totally non-singular B is totally non¬ 
singular as well: Assume for a contradiction that B is not totally non-singular, i.e. that B 
has a vanishing entry. By permuting rows we can assume A and B to be of the form 



for some matrices F G M{l,q), D G M{q,q) and E £ M(q, 1). But then 

A*B = F*0 + D*E = D*E = 0. 

In particular this implies that the q x q submatrix D of A is singular, contradicting the 
fact that A is totally non-singular by assumption. 

Induction step. Assume the claim holds for an n > q and let A G M(n -|- l,g) be 
totally non-singular. Note that for each i G {0,..., n}, the n x q matrix A^j is totally non¬ 
singular since A is totally non-singular. Thus, by the induction hypothesis, we can find for 
each i G {0,..., n} a totally non-singular matrix Ci G M{n, n — q) such that A* Ci = 0. 
For f G {0,..., n},j £ {0,..., n — q} let C{i, j) £ M{n + 1, n 1 — q) be the matrix with 
C{i, j)’l = Ci and 0 else. Then, for alH G {0,..., n}, j G {0,..., n — g}, C{i,j)’l is totally 
non-singular, C{i,j)[j] = 0 and A*C{i,j) = 0 by construction. 

Step 1. First, for each i £ {0,...,n}, we deform C{i,0) into a matrix C{i,0) £ 
M(n -£ l,n -\-1 — q) with the following properties: 

1. A*C'(z,0) = 0, 

2. C{i,0)'^ is totally non-singular, 

3. All {n-\- 1 — q) x {n-\- 1 — q) minors of C{i, 0)^i are nonzero. 
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Let i G {0,..., n}. For a := (/cq, ..., /c„-q) G S := {(/q, • • •, In-c^ : 0 < < • • • < 

In-q <n — l} define the projection : M{n + l,n + 1 — q) ^ M{n + 1 — q^n + 1 — q) hy 
Pa{X){j} ;= {X^i){kj} for all X G M(n + 1, n + 1 — g), j G {0,..., n — g}. Now let a G S, 
and set := Pcr{C{i,0)). By permuting rows we can assume A and C{i,0) to be of the 
form 

C(i.0)=(®;) (7) 


for some matrices F G M{n + 1 — q, q), D £ M{q, q) and F^^ G M{q, n + 1 — q). 


Next, we show that there is a vector u„ = 


where G G 


such that A*Ucr = 0 and det(i?o- + Paiu^CQ)) = det(-Eo- + VctCq) 7^ 0: Since Ci is totally 
non-singular, Fa- has rank n — q. Thus we can find a vector G such that 

det(-Eo- + ^^0-^0) 7^ 0 Finally, we just have to ensure that A*Ua = 0. Since D is 
totally non-singular there is a vector Wa G C? such that D*Wa = —F*Va and this gives 
V 

= _p*r;o- + D*Wa = 0. Repeating this construction, we can find a collection of 


A 


Wr 


vectors {tto-jo-es ^ such that for all a G S we have A*Ua = 0 and det [Pa{C{i, 0) -|- 

'^o-Gq)) 7 ^ 0 . 

Next, for distinct cri,a 2 G S, define the mapping K{X) := C{i,0) + u^^eQ -|- 
A G C and note that by construction A*K{\) = 0 for all A G C. Note that by construction 
K{X)’^ = Ci is totally non-singular for all A. Furthermore, the (n -|- 1 — g) x (n -|- 1 — g) 
minors det(Po-i (R"(A))) and det(Po-2 considered as polynomials in A. The 

polynomial equations det(Po-i (-f^(A))) = 0 and det(Po-2(iL(A))) = 0 are non-trivial: For 
A = 0 we have det(Po-i (R"(0))) = det(Po-i 0) + ^o-iCg)) / 0 by construction. For A 
large one can consider jUa^eg as a small perturbation to Uo-jCg. Thus, using linearity of 
the determinant in the 0-th column, we conclude that 


det(Po-2(iL(A))) = det(PCT2(C'(bO) + Ua^el -F A^aaeo)) 

= A •det(Po-2(C'(i,0) + Ua^e^) + jPa^iuaiel)) / 0 

for large enough A by the continuity of the determinant and the fact that det(Po-2(C*(i, 0) -|- 
^0-260)) 7^ 0 construction. A non-trivial polynomial equation in one variable just has a 
finite set of solutions and hence the set 


{A G C : det(P<,i(iL(A))) = 0 V det(P^2(^(A))) = 0} 

= {A G C : det(P<,i(A:(A))) = 0} U {A G C : det(P^2(^(A))) = 0} 

is finite. In particular there is an ao-2 G C such that det(P(7^(Ar(ao-2))) 7^ 0 and 
det[Pa 2 {K{ 0 , 0 - 2 ))) 7^ 0 Applying the same argument to L{X) := C{i,0) -|- UoiOq -|- 


® The direct sum decomposition of Uo is with respect to the decomosition given by Equation (7), i.e. 
A*Uo = F*Vo + D*Wa. 

Note that i?o-[0] = 0 by construction of C{i, 0). 

In fact this holds for almost all 0^2 G C. 
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aa 2 UcT 2 ^Q + Atto-geQ, A G C, where 1 T 3 G S is distinct from (Ti,(T 2 , yields an Oo-g G C 
such that det[Pai{L{aa 3 ))) 7 ^ 0 , det(Po- 2 (-L(ao- 3 ))) / 0 and det(Po- 3 (L(a(T 3 ))) 7 ^ 0 
Finally, since |S| is finite, we can inductively apply the argument to obtain a matrix 
C(i, 0) = C{i, 0) + Mo-iCq + Xlo-eS CLaUaSQ with the desired properties. 

Step 2. Secondly, we construct for each i G {0,..., n} a matrix Di G M(n + 1, q) with 
the following properties: 

1 . A*Di = 0 . 

2 . {Di)^i is totally non-singular. 

Let i G {0,...,n}. Let Di{Xi,..., Xn-q) '■= C{i,0) + j) where Xj G C, 

j G {1,..., n — q}, and note that by construction we have A*Di(Xi ,..., Xn-q) = 0 for 
all Ai,..., Xn-q G C. By choosing (Ai,..., Xn-q) appropriately one can make sure that 
{Di{Xi,..., Xn-q)),i is totally non-singular: First let G{X) := C{i,0) + XC{i,l), A G C. 
Just like in Step 1, the minors of G'(A)’^ and G(X)’j together with the (n + l — q) x (n+l — q) 
minors of G(X)^i yield a finite set of polynomial equations in A. All of these polynomial 
equations are non-trivial: For A = 0 none of the minors of G'(O)’^ = C* and none of the 
(n + 1 — q) X (n + 1 — q) minors of G(0)^i = G(z,0)^i vanish by construction of G(z,0). 
For large A one can consider ^(7(i,0) as a small perturbation to G(i, 1). Hence, for large 
enough A, none of the minors of jG(A)’| vanishes by the fact that C(i, 1)’| = Gi is totally 
non-singular by construction and the continuity of the minors. Thus, just like in Step 1, we 
conclude that there are just finitely many values of A for which any of these polynomials 
vanishes. In particular there is an Ai G C such that both G(Ai)’^ and G(Ai)’| are totally 
non-singular and all {n + l — q) x {n + l — q) minors of G(Ai)^j are nonzero. Applying the 
same argument to H{X) := C'(i,0) + XiG{i,l) + AC'(i, 2 ), A G C, yields an A 2 G C such 
that H (A 2 )’^, H (A 2 )’j and H (A 2 )’^ are totally non-singular and all {n + 1 — q) x {n+1 — q) 
minors of iL(A 2 ),i are nonzero. Choosing the values for Xj, j G {1,... ,n — q}, inductively 
in this fashion finally yields a matrix Di with the desired properties. 

Step 3. To complete the induction step we choose by a similar argument as in Step 
1 and Step 2 before parameters 7 ^- G C, j G {1,... ,n}, in B := Dq + such 

that Bj is totally non-singular for each i G {0, ...,n}, i.e. such that B is totally non¬ 
singular: First define 1(A) := Dq + XDi, A G C. Clearly /(0)^o = (-Oo),o is totally 
non-singular by construction of Dq. Furthermore, for large A , jDq can be considered 
as a small perturbation to Di. Thus, for A large enough, jI{X)^i is totally non-singular 
by construction of Di and the continuity of the minors. Hence, all the minors of L(A)^o 
and /(A )7 yield non-trivial polynomial equations in A and therefore there are just finitely 
many values for A for which any of these minors vanishes. In particular there is a 71 G C 
such that both /( 7 i),o ^-re totally non-singular. Applying the same argument 

to J(A) := Dq J- 7 iiJ 2 + AiJ 2 yields a 72 G C such that 7 ( 72 ),O) «^( 72 ),i and 7 ( 72),2 are 
totally non-singular. Continuing to choose the 77 i G {1,..., n}, inductively in this fashion 
then yields a totally non-singular matrix B with A*B = 0. □ 
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Also in this case we obtain a finite set of non-trivial polynomial equations in A and thus the argument 
given before can be applied to find 00 - 3 . 
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Lemma V.2. Let g G {1,..., n — 1} and let A £ M{n,q,M.) be totally non-singular. Then, 
there exists a totally non-singular matrix B G M(n, n — g,M) sueh that A^B = 0. 

Proof. The arguments given in the proof of Lemma V.l also apply to real numbers. □ 

For A: G {1,..., 2n —3}, define the inclusion in the k-th. antidiagonal H{n) 

by 


^ ( Vj A j + l = k,j <l 

{t'k{v))ji v*i A j + l = k,l < j 

\ 0 else 

where 

, ,. f \k/2\ A k <n — 1 

^(".'=) = (f„_l_t/21 iffe>„-l 

is the length of the upper half of the A:-th antidiagonal. By expanding in the generalised 
Gell-Mann orthonormal basis of H{n), it is easily seen that the inclusion of real vectors in 
the same antidiagonal preserves the standard inner product, i.e. for k G {1,... ,2n — 3} 
we have 

tT:{^Lk{v)ik{w)) = {v,w), yv,w G (8) 

Furthermore, the inclusion of an imaginary and a real vector in the same antidiagonal 
yields Hilbert-Schmidt orthogonal matrices, i.e. for k G {1,..., 2n — 3} we have 

tr(rfc(u)tfc(ra)) = 0, \/v,w G (9) 

and finally that inclusions of vectors in different antidiagonals also yield Hilbert-Schmidt 
orthogonal matrices, i.e. for k,j G {1,..., 2n — 3} with k ^ j we have 

tr{ik{v)tj{w)) = 0, Vu G u; G (10) 

The following theorem is the main result of the present paper. 

Theorem V.3. Let Gq be a basis of Cf . Furthermore, for k £ {2r— 1,... ,2{n — r) — T}, 
let Ak,A'i^ £ M('y{n,k),r,R) be totally non-singular and define the tuple 

Gk := {Lk{Ak[0]),LkiiA'j^[0]),ikiAk[l]), Lk{iA'k[l]),.. .,Lk{Ak[r - 1]), Lk{iA',^[r - 1])). 

Then G := GqU G 2 r-i U G 2 r U • • • U G 2 (n-r)-i is r-eomplete and |G| = 4r(n — r) -|- n — 2r. 


C" was defined in Equation (4). 
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Proof. The idea of the proof is to use Lemma V.2 to determine a basis of the null space 
of Mg such that the construction of [19] can be applied. We do this in the first step of 
the proof. In the second step of the proof we use the construction of [19] to show that Mq 
indeed is r-complete. 

Step 1. First, by Lemma V.2, there are totally non-singular G 

M{j{n, k),j{n, /c) — r, M), k e {2r -|- 1,..., 2(n — r) — 3}, such that 


— 0 , 

= 0 . 


( 11 ) 


Now let 


Gi ■■= (tfc(Bfc[0]),ifc(i5fc[0 ])>---- r- 1]), ik{iB'ffa{n,k) - r- 1])), 

k G {2r -1- 1,..., 2(n — r) — 3} 

and let Gj: = (0) for k G {2r — 1, 2r, 2(n —r) — 2, 2(n —r) — 1}. In the remainder of this first 
step we prove that U • • • U is a basis of Ker(MG'): For /c G {1,..., 2n — 3}, 

let Qk '■= {X G H{n) : X{j) = 0 Vj / /c A Xu = 0 for 2i = k} be the subspace of 
Hermitian matrices with vanishing diagonal and non-vanishing entries only in the A:-th 
antidiagonal. By Equation (10), H{n) can be decomposed into the following mutually 
orthogonal subspaces: 


Hin) — Cf © Q2r-1 © • • • © Q2(n-r)-l- (12) 

Note that Span(GA: U G^) C Qj. for all k G {2r — 1,..., 2(n — r) — 1}. Hence, by the 
decomposition (12), to show that G^r-i U • • • U is a basis of Ker(MG) it suffices 

to prove that for k G {2r — l,...,2(n — r) — 1} the matrices G^ U Gk span the subspace 
Qk and that Span(G^) C Ker(MG'). First observe that indeed Span(G^) C Ker(MG') for 
every k G {2r + 1,..., 2(n — r) — 3}: Note that for every k G {2r + 1,..., 2(n — r) — 3}, 

iT{ik{Ak[l])ik{Bk[j])) = {A-k[l],Bk[j]) = {AlBk)ij = 0, 
tr{ik{iA^,[l])ik{iBl[j])) = {A'k[l],Bi[j]) = {{A^k)'Bi)ij = 0, (13) 

V/ G {0,..., r - 1}, j G {0,..., 7 (n, k) - r - 1}, 

by equations (8) and (11). Furthermore, 

ti {ik{iA'k[l])ik{Bk[j])) =0, 

ii{ik{Ak[l])ik{iB'k[j])) (14) 

V/ G {0,..., r - 1}, j G {0,..., 7 (n, k) - r - 1}, 

by Equation (9). I.e. Span(G^) is orthogonal to Span(Gfc) and thus Span(G^) C Ker(MG'). 

To conclude the first step, we prove that G^ U Gk spans the subspace Qk for every 
k G {2r — 1,..., 2(n — r) — 1}: Let k G {2r — 1,..., 2(n — r) — 1}. Since Ak is totally 
non-singular, the columns of Ak are linearly independent and the same argument applies 
to A'j.. Hence, by the equations (8) and (9), Gk is a tuple of linearly independent Hermitian 
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matrices. The same argument applies to k G {2r+l,..., 2(n—r)—3}. But we have seen 
that Span(Gfc) is orthogonal to Span(G^) for k G { 2 r + 1 ,..., 2(n — r) — 3}. Furthermore, 
for /c G {2r - 1,..., 2(n - r) - 1}, |G^| + \Gk\ = k) -r) + 2r = 2^{n, k) = dimQfc 

and thus Gj: U G^ indeed spans Q^. 

Finally, observe that 


2(n—r) — l 2r—2 2{n—2r)+l 

|G|=dimG”+ ^ |Gi| = ^ 27 ( 71 , i) + n + ^ 2r 

i=2r—l 2=1 2=1 

= (2r)^ — 2(2r) + n + 2r{2{n — 2r) + 1 ) 

= 4r(n — r) + n — 2r. 


Step 2. In the second step, we essentially reproduce the construction of [19] and some 
ideas of [26]. We show in the following that every nonzero matrix X G Ker(MG') has at least 
r + 1 positive and r + 1 negative eigenvalues and this concludes the proof by Proposition 
II.2. 

Let X G Ker(MG') be arbitrary. By the interlaced eigenvalue Theorem (Theorem 4.3.15 
of [28]) it suffices to prove that there is an 2 (r + 1 ) x 2 (r + 1 ) principal submatrix of 
X with r + 1 positive and r + 1 negative eigenvalues. We conclude the proof by finding 
such a submatrix: There is a smallest number k G {2r + 1,..., 2(n — r) — 3} such that X 
has non-vanishing entries in the fc-th antidiagonal. First note that either the real or the 
imaginary part of the A:-th antidiagonal does not vanish. Let us consider the case where 
the real part does not vanish, the other case can be shown analogously. The real part of 
the /c-th antidiagonal of Ker(MG) is spanned by the 'y{n,k) — r real matrices of G^, i.e. 
each X G Ker(MG') is a linear combination of the j{n,k) — r real matrices of G^. But 
then there have to be at least 2(r -|- 1) non-vanishing entries in the fc-th antidiagonal of X 
because otherwise there would be a vanishing ( 7 ( 77 , k) — r) x (j(n, k) — r) minor of Bk and 
this contradicts the fact that Bf^ is totally non-singular (For more details see Lemma 9 of 
[26].). I.e. there is a 2(r -|- 1) x 2(r -|- 1) principal submatrix of X of the from: 


/O 

0 

0 .. 

0 

0 

Xl 

\ 

0 

0 

0 .. 

0 

X2 

tl 


0 

0 

0 .. 

. X3 

yi 

yl 


0 

0 

X3 . . 

0 

y2r-2 

t2r- 

1 

0 

X2 

yl ■■ 

• ylr- 

2 0 

Vlr 


\xi 

yl 

yl ■■ 

■ ylr- 

■1 ylr 

0 

/ 

G C, 

3 e 

{ 1 ,.. 

.,r },7 

G {1,... 

,, 2 (r 

+ 


,XiGC\{0}, 7 G {1,..., r-|-1}, 


(15) 


where yl 

Finally, we show by induction that a matrix of this form has at least r + 1 positive 
and r -|- 1 negative eigenvalues: The claim clearly holds for r = 0. Now assume the claim 
holds for r G Nq. Let T be a 2(r -|- 2) x 2(r -|- 2) matrix that is of the form illustrated in 
Equation (15). Then, one can obtain a principal 2(r -|- 1) x 2(r -|- 1) submatrix Y' of Y 
that is of the same form by e.g. deleting the first and last row as well as the first and last 
column of Y. Thus, by the induction hypothesis and the interlaced eigenvalue Theorem 
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(Theorem 4.3.15 of [28]), Y has at least r + 1 positive and r + 1 negative eigenvalues. A 
straightforward calculation shows that det(y) • det(y^) < 0. Since the determinant of a 
matrix is the product of its eigenvalues, the claim follows from det(y) • det(y') <0. □ 

In the following the r = 1 case is of particular interest because Theorem III.l is obtained 
from this case by choosing the totally non-singular matrices appropriately. 

Corollary V.4. Let Gq := (coCq, ..., en-ie*_;^). Furthermore let Wk,Vk G k G 

{1,..., 2n — 3}, be sueh that every entry of Vk and every entry of Wk is nonzero. Then 
G := Go U U ••• U {i2n-2.{v2n-z), i^2n-z{iw2n-z)) is l-complete and |G| = 

5n-6. 


Proof. First, note that Gq is a basis of Gf. Furthermore as by assumption all en¬ 
tries of matrices Wk,Vk G ~ k G {l,...,2n — 3}, are nonzero, 

we conclude that all their minors are nonzero^^. Consequently the matrices Wk,Vk G 
M( 7 (n,/c), 1,M), k G {1,..., 2n — 3}, are totally non-singular. Hence Go and Gk '.= 
(ii(tJi), ii(rai)), k G {1,..., 2n — 3} fulfil the conditions of Theorem V.3 for r = 1 and 
thus G = Go U Gi U ... U G 2 n -3 is l-complete. □ 

Example For i G {1, ... ,2n — 3}, we can choose Wi = Vi = ^/2e, where e := (1,..., 1) G 
is the vector with a one in every component. Altogether this yields 2(2n — 3)+n = 
5n — 6 Hermitian operators for G. For n = 4 these are 


/I 

0 

0 

0\ 


0 

0 

0\ 



0 

0 

0\ 


0 

0 

0\ 



1 

0 

0\ 

( ° 


0 

0\ 


0 

1 

0\ 


(° 

0 


0\ 

0 

0 

0 


0 

1 

0 

0 


0 

0 

0 

° 

0 

0 

0 

0 


1 

0 

0 


— i 

0 

0 


0 

0 

0 

0 


0 

0 

0 

° 

0 

0 

0 

0 ’ 

0 

0 

0 

0 

’ 

0 

0 


0 ’ 

0 

0 

0 

0 

’ 

0 

0 

0 

0 ’ 

0 

0 

0 

0 ’ 

1 

0 

0 

0 

’ 

— i 

0 

0 

0 

\o 

0 

0 

0/ 

^0 

0 

0 



^0 

0 

0 

0/ 

^0 

0 

0 



^0 

0 

0 

0/ 


0 

0 

0/ 

^0 

0 

0 

0^ 



0 

0 

0/ 



0 

0 

1\ , 

( 0 

0 

0 

A 1 


0 

0 

0\ , 


0 

0 

0\ , 


0 

0 

0\ 



0 

0 

0 

0 

1 


0 

0 

i 

° 

0 

0 

0 

1 

0 

0 

0 

M 

0 

0 

0 

0 

1 1 

0 

0 

0 

0 

1 

0 

0 ’ 

0 

— i 

0 

O ’ 

0 

0 

0 

0 ’ 

0 

0 

0 

0 ’ 

0 

0 

0 

1 

’ 

0 

0 

0 


0 

0 

0/ 


0 

0 

0/ 

^0 

1 

0 

0 

^0 

— i 

0 

0/ 

^0 

0 

1 

0^ 


^0 

0 

— i 


Finally, let us give a proof of Theorem III.3. 

Proof. For k G {2r — l,...,2(n — r) — 1}, define Ak,A'^ G M( 7 (n,/c), r, M) by setting 
{Ak)ji = {A'^)ji = for all j G {0,..., 7 (n, A:) — 1}, ^ G {0,..., r — 1}. Observe that 
both Ak and A'^ can be considered as the first r columns of a 7 (n, k) x 7 (n, h) Vandermonde 
matrix and since Xj ^ xi for all j, Z G { 1 ,..., r} with j ^ I and xi ^ 0 for all Z G { 1 ,..., r} 
they are thus totally non-singular. Applying Theorem V.3 to Ak,A'^ then concludes the 
proof. □ 


B. Proof of Theorem III.l 

Let us now give a proof of Theorem III.l. 
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The minors of a m x 1 matrix are simply the entries of the matrix. 
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Proof. Define G H{n), A: G {1,..., 2n — 3}, by 

1-1 


■— <^j+Z,fc COS 
{Yk)ji := i5j+z,fcsin 


-TT 


2n 

3-1 ^ 

2n 

j,l e {0,...,n- 1}. 


-TT 


Next observe two things: 


1. The matrices {Xi, Yi,..., X 2 n- 3 , Y 2 n- 3 } T YA(n) are linearly independent by equa¬ 
tions (9) and (10). 


2. Since 0 < < f for j, / G {0,..., n — 1}, j > /, we hnd {Xk)ji / 0 and (Yfc)j 7 / 0 

for j + I = k, j > 1. 

Let Uk,Wk G /c G {1,... ,2n-3}, be such that ifc(rtfc) = - 5fc/2,|-fc/2]e|-fc/2]ep^/2-|, 

I'kiiwk) = Yf^. and note that both and have no vanishing entry. Thus, by Corollary 
V.4, G := (eoco,... ,e„_ie*_i) U (Xi, Yl, ..., X2n-3,Y2n-3) is 1-complete. 

Let G := {eoe^,... ,en-ie^_i,vivl,vivl,... ,V2n-3V2n-3X2n-3V2n-3)- To conclude the 
proof, we show that Span(G) = Span(G). First note that for k G {1,..., 2n — 3} 


2 n—3 

J = 1 
2n—3 

Vkvl = Y, xiiXj - Y,) + eoeS + 
i=i 


and thus Span(G) C Span(G). In order to show that Span(G) C Span(G), consider the 
matrix 


/ Xi 

^2 

Xi 

rf,3 

Xl 

^2n 3 
. X^ 

X2 


™3 

2 n—3 

^2 

X 2 

• ^2 


2 

3 

2n—3 

X3 

^3 

3^3 •• 

• 2:3 

\X2n-3 

^L-3 

r^3 

X2n-3 ■ ■ 

2n-3 
■ X2n-3- 


The matrix T is a Vandermonde matrix and thus invertible if Xi ^ xj for i ^ j. Hence we 
hnd 

^ 2n—3 

= 2 ^ - 2eoe5 - 2x1^-"^en-iel_^), 

i=i 

2n—3 

i=i 
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Different from the rest of the present paper, the indices we use to label T begin with 1, not with G. 
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and this shows that Span(G) C Span(G). □ 

Remark Note that there are many possible choices for the phases of the Vi. The only 
constraint is that j(p ^ ^ for all j G {1,..., n — 1}, /c G Z. 

C. Proof of Theorem IV. 1 and Proposition IV.2 

For Xy G S^, E G H{n), e > 0 and a measurement M define the set 

F,{Xr, E, M) := {y G 5- : \\M{Y) - bh < e}, (16) 

where b = M{Xr + E). 

Lemma V.5. (Stability.) Let M be an r-complete measurement and let e > 0. Then, 
there exists a constant Cm > 0 independent of e such that for all Xj. G Sf, and E H (n) 
with ||M (£')||2 < e we have 

Y G F,{Xr,E,M) => ||y-< Cme- 

Proof. Denote by vr : H{n) —s- Range(M*) the orthogonal projection on the subspace 
Range(M*) C H{n) and by tt-*- ; H{n) —)• Ker(M) the orthogonal projection on the 
subspace Ker(M) C H{n). Furthermore, let Y' := 'K-^{Xr) + 7r(y) and let amin be the 
smallest singular value of M Then, we find 

||X, - y'li = ||7r(v, - y )||2 < —\\M{Y - V ,)||2 

^min 

< —{\\M{Xr) - bh + \\M{Y) - bh) < —(||M(E)||2 + e) 

^min ^min 

- 

^ min 

From the spectral variation bound for Hermitian matrices (Theorem III.2.8 of [29]) we 
conclude that 

||Eig(X,)-Eig(y')||2 = 

< —e 

^min 

But this implies that I ART') I < —for i G |r + 1,... ,n|. 

^ ^ '^min ^ -* 

Next, note that 

K := — max Xn-r(Z) 

ZeKer(M),||Z||2=l 
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We assume M to have full rank. 
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exists by compactness of {Z G Ker(M) ; \\Z\\‘2 = 1} and continuity of \n-r- Furthermore, 
by Proposition II.2, every nonzero Z G Ker(M) has at least r + 1 negative eigenvalues and 
hence we conclude that k > 0. 

There exists Z G Ker(M) with \\Z\\2 = 1 and a > 0 such that Y = Y' + aZ . Since 
y > 0 we conclude from Weyl’s inequality (Theorem III.2.1 of [29]) that 

0 < \n{Y' + CtZ) < Ar+l(yO + Ol\n-r{Z) < -C — UK. 

^min 

and hence we find 

2 , ^ 
a < -e. (18) 

Finally, combining equations (17) and (18), we conclude that 


|F - X,||2 = ||y' + aZ- X,||2 < ||y' - X,||2 + lla^lb 


< 


+ 


(Tri 


Ka„ 


e. 


Choosing Cm = (1 + -) then proves the claim. 


□ 


Remark Since k just depends on Ker(M), it is independent of the choice of basis for 
Range(M*). Thus, since it is always possible to choose an orthonormal basis of Range(M*), 
the constant Cm is mainly determined by k. 


The proof of Theorem IV. 1 is an immediate consequence of this lemma. 

Remark Let M be a measurement that is not r-complete. Then there exist Z^ G and 
Z G with Zr ^ Z such that M{Zr — Z) = 0 and we find Z G F^{Zr, E, M) for all e > 0 
and E G H{n) with ||M(F1)||2 < e. Thus, if 1 G Range(M*), the r-complete property is 
necessary to enable the recovery of every X^. G via the optimization problem (5). 

Finally let us give the proof of Proposition IV.2. 

Proof. From Theorem IV.1 we obtain the bound [[V — xx *||2 < Cm^ and the proof of 
Lemma IV yields the bound \/Y17=2 — C'm^. From this we find 

||xx* — XX*||2 < ||y — XX*||2 + ||v — XX*||2 < 2Cm€- 


Note that aZ = n^iY) - Ti^{Xr) G Ker(M). 
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Finally let (p G [0, 27r) be such that (x, e*‘^x) is positive. Then, 


X 

|2 -P 

l2 T 

X 

|2 
l2 ■ 

- 2Re ((x, e^'^x))) x 

||2 

Il2 


X 

|2 -P 

l2 T 

X 

|2 
l2 ■ 

1 

ro 

IT 

to to 



X 

|2 -P 

l2 T 

X 

|2 
l2 ■ 

1 

to 

IT 

to to 

+ 

^lli 

+ 2|(x,x)|) 

X 

li + 

X 

li) 

to 

1 

IT 

to 



X ; 

^ + 11 

X [ 

1 

2 

2 (x, x) p + 2 X 2 X 

||2 

Il2 

- 2|(x,x)p 


< 2 (||x||| + ||x||^ - 2|(x,x)|^) 
= 2||xx* — XX *||2 

< 2{2CMef- 
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